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Abstract 



We present the full one-loop calculation for gg — > ZZ in the Minimal Su- 
persymmetric Model (MSSM) including nonresonant contributions from the 
squark loop diagrams and provide analytical expressions for the helicity am- 
plitudes. The one-loop process gg — > ZZ via quark loops has been calculated 
in the Standard Model. In supersymmetric models, additional contributions 
arise from squark loops. In some regions of the MSSM parameter space, the 
top and bottom squark loops can make important contributions to the dia- 
grams involving Higgs bosons. The heavy Higgs scalar (H) might be detected 
at the Large Hadron Collider via gg — > H — > ZZ for tan j3 ^ 5. 
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I. INTRODUCTION 



One of the main purposes for constructing a hadron supercollider is to discover the 
mechanism responsible for the spontaneous breaking of the electroweak symmetry. A large 
amount of attention has been devoted recently to the ability of detecting a Higgs boson at 
the CERN Large Hadron Collider (LHC). Gluon fusion is the main mechanism of producing 
Higgs bosons at the LHC, and if a CP-even Higgs scalar is sufficiently heavy it can decay 
into Z boson pairs. In the Standard Model, this is a useful detection channel for the Higgs 
boson with a mass M HsM > 2M Z since the HsmZZ coupling grows with Higgs boson mass. 

The physical Higgs bosons of the minimal supersymmetric extension to the Standard 
Model (MSSM) is comprised of two CP-even states, a lighter h and a heavier H, one CP-odd 
state, A, and two charged Higgs bosons H ± . The Higgs potential is constrained by super- 
symmetry such that all tree- level Higgs boson masses and couplings are determined by just 
two independent parameters, commonly chosen to be the mass of the CP-odd pseudoscalar 
(mA) and the ratio of vacuum expectation values (VEVs) of Higgs fields (tan/3 = vijv\). 
The CP-even states can couple to the ZZ final state at the tree level, but only the heavier 
one (H) can be massive enough to decay into on-shell Z bosons. For m# ~ rriA ~ 2m 4 , 
discovery of the MSSM heavy Higgs scalar in the ZZ channel could be problematic due to 
the fact that the Higgs couplings to the gauge bosons arise from the D-term contributions to 
the scalar potential. As the H becomes heavy, its coupling with the gauge bosons becomes 
smaller. Previous studies have found that the production rate for gg — > H — > ZZ is large 
enough to make this channel a viable one for Higgs detection provided tan/5 = t^/fi is 
relatively small 5) and that the Higgs boson is not too heavy (m# < 350 GeV) 
The impact of supersymmetric decay modes has also been investigated ||. 

The production cross section of gg — > (0 = H,h,A, and Hsm) has been studied 
in increasing detail in the last few years, with much attention devoted to the Standard 
Model Higgs boson or to the lighter MSSM Higgs scalar, h. The impact of squark loops 
in the production process has been included, and the impact of squark mixing has been 
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investigated |J. The QCD corrections to the production process have been calculated and 
shown to give a large enhancement to the cross section 0. 

It is well known that the dominant background to the channel gg — > Hsm — * ZZ 
in the Standard Model is qq — > ZZ, and the one-loop continuum production of Z pairs via 
quark loops, gg — > ZZ. In the MSSM the existence of the scalar partners to the quarks 
gives a new contribution to the continuum production that can interfere with the quark loops 
and possibly affect the overall level of background. If the masses of the squarks are much 
larger than the electroweak scale, their contributions to the signal should be small as their 
mass arises from soft terms rather than from the Higgs mechanism. For squarks with masses 
comparable to the electroweak scale, significant contributions to the Higgs production signal 
might be expected. It has been shown |5|j9|JT0|] that the squark loops can have a significant 
impact on the production cross section if the squarks are fairly light. Therefore, a natural 
step is to also investigate the impact of the squark loops on the background. 

The Feynman diagrams for the background coming from squark loops to the Higgs signal 
in the channel gg — > H —>■ ZZ is shown in Fig. 1. Since only the third generation gives a 
sizable contribution to the resonant Higgs contribution gg — > H — > ZZ, one might expect 
an important contribution from the continuum background for which all three generation 
squarks are expected to contribute. The central issue concerning the detectability of the 
Higgs boson is whether the squark loop contribution can be comparable to the quark loop 
contribution for energies near the Higgs mass. The squark loop contribution is maximal near 
the production threshold (y/s = 2m^), so we are most interested in the size of the squark 
loop background for 2m 9 - m Mr. 

Another potential application of the process under consideration is a measurement of the 
Hgg coupling. The production mechanism involves a coupling to the Higgs boson (a Yukawa 
coupling) that implies that only the third generation matter make a sizeable contribution 
to the rate. Each generation contributes equally to the background (assuming equal masses 
for squarks in different generations), so potentially there is an enhancement as there is in 
the case of quark loops. A large background from squark loops could conceivably alter 
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the regions of parameter space (the Ma — tan (3 plane) for which the heavy Higgs can be 
discovered. 

The gluon coupling to the Higgs boson is mediated by triangle graphs. The calculation 
of the background involves box graphs and is significantly more complicated than that of 
the Higgs production cross section. The one-loop integrals, which are expressible in terms of 
Spence functions (dilogarithms) ||TTl , |l2l , were evaluated numerically using a FORTRAN code 
P~3|1 . We have not included squark mixing in this first survey, and leave the more general 



case to a future paper ||14|| . However we do not expect squark mixing to have a significant 



effect on the cross sections. While the cross section for the resonant Higgs graphs can be 
significantly increased by an enhanced Htiti coupling [|IIJ], squark mixing results only in 
mixing angle factors in the couplings of squarks to the Z bosons and splits the degeneracy 
of the squark masses for the nonresonant graphs. One expects an enhanced contribution 
only from the loops involving t\ squarks which might have a suppressed mass with significant 
mixing. If the lightest squark is less than half the mass of the Higgs boson, then the Higgs 
boson can decay into squark pairs making the Higgs broader and more difficult to detect in 
the ZZ final state. 
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Fig. 1(a). The bubble diagrams. 
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Fig.l(b). The triangle diagrams. 
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Fig.l(c). The box diagrams. 
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Fig.l(d). The diagrams involving the Higgs bosons. 
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II. THE SQUARK LOOP CONTRIBUTION TO THE BACKGROUND 



In this section we compare the size of the squark loop background with that arising from 
the quark loops. In the Standard Model, only quark loops contribute to the production of ZZ 
from gluon fusion gg — > ZZ. In the MSSM there are additional contributions from squark 
loops, which typically are thought to be suppressed by large squark mass. In this section 
we present the cross sections expected at the LHC with yfs = 14 TeV. In our numerical 
computations we use the CTEQ3 parton distribution functions [[HJ with A 4 = 0.177 GeV. 
We take a" 1 = 128, sin 2 6 W = 0.2319, a s {M z ) = 0.124 and m t = 175 GeV. 

In Figure 2 we plot the invariant mass distribution da/dMzz{pp — ► ZZ + X) from 
gg — > ZZ via the squark loops in the MSSM, assuming their masses are all degenerate at 
values of 100 GeV, 300 GeV and 500 GeV. The lightest squark considered is excluded by the 



Tevatron data HIE] , but we include it for comparison purposes. There is a rapid rise in the 
cross section until the threshold of real squark pair production after which the cross section 
falls again. Since the squark loop cross section peaks at y/s ~ 2m,, the largest contribution 
to the background from the squarks comes at subprocess {gg — > ZZ) energies slightly below 
the squark production threshold. The contribution to gg — > ZZ from graphs with squark 
loops not involving the Higgs boson is separately gauge invariant, and the cross section falls 
off rapidly as y/s increases. 

In Figure 3 we show the squark loop cross section for the polarization states LL, TT and 
TL + LT of the ZZ final state. The contribution from squark loops is always smaller than 
that from quark loops even for squarks as light as 100 GeV. The transversely polarized case 
dominates over the other modes for the squark loops just as it does for the contributions 
from the quark loops. 
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Fig. 2. da/dMzz{pp ZZ + X), from gg — ► ZZ, via squark loops with 12 squarks, 
without Higgs bosons, for nig = 100 GeV, 300 GeV and 500 GeV. Also shown arc 
the contributions from the box diagrams of quark loops. These contributions are 
independent of tan/3 for fixed squark mass rriq. 
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Fig. 3. da/dMzz{pp —> ZZ + X), from gg — > ZZ, via squark loops with 12 squarks, 
without Higgs bosons, in the (a) LL, (b) TT and (c) TL + LT states of the ZZ for 
trig = 100 GeV, 300 GeV and 500 GeV. 
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III. THE SUPERSYMMETRIC LIMIT 



Some contributions to the one-loop amplitudes for squark loops given in the Appendix 
are closely related to parts of the helicity amplitudes for the quark loops. These types of 
super symmetric relationships Jl7|| were observed in the explicitly computed weak interaction 
process Z — > 37 or equivalently 77 — > 7Z fllSfl , and for 77 — > 77 fl9fl . One can also compare 



the contributions from fermion loops in the process 77 — > HH |20] to the contribution that 
was obtained later for the gauge boson loops [[Hfl. However the most numerically significant 
contributions to the cross section is often those contributions that are not given by the 
super symmetric relationship. In fact the asymptotic behavior (y/s — > 00) is governed by the 
spin of the particles exchanged in the t- and w-channels. For example, one has the dominance 
of the W boson loops over fermion loops in the process 77 — > ZZ |22[ . 



14 TeV, gg -> ZZ 



>10" 
a 



10" 



x 
+ 

ISI 



t 10" 

a, 
a, 



30 

\ 

b 

T3 



10 



10 



-12 



1 1 1 1 1 1 1 1 1 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 
(c) TL+LT 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



200 400 600 800 200 400 600 800 200 400 600 800 1000 
M zz (GeV) M zz (GeV) M zz (GeV) 



Fig. 4. Differential cross section of pp — > ZZ + X from gg — > ZZ via top quark and 
top squark loops without Higgs bosons for m; = =175 GeV. 

For the process gg — > ZZ discussed in this paper, the SUSY relationship exists between 
the squark loop amplitudes and the parts of the quark loop amplitudes arising from the 
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vector coupling of the Z boson to quarks. This relationship provides another consistency 
check on our analytic results in addition to the usual check of gauge invariance. In Fig. 4 
the contributions from the top quark loops is compared to the contributions from II and tR 
loops for masses set equal to m t = m~ t = 175 GeV. 

IV. INTERFERENCE EFFECTS 

In Fig. 5 we present the contributions to gg — > ZZ including both the resonant Higgs 
graphs and the nonresonant graphs for the quark loops alone and for the squark loops alone 
with tan (3 = 2 and 10. The squark mass is 200 GeV, and M A is 300 GeV for which M H ~ 309 
GeV. For tan (5 > 5 the differential cross section appears the same apart from the size of the 
Higgs boson resonance. 



vs" = 14 TeV, tan0 = 2, M A = 300 GeV 
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Fig. 5. Differential cross section of the Higgs signal in pp collisions from (a) quark 
loops and (b) squark loops. The dashed lines are the contributions from the Higgs 
graphs, the dotted lines are the contribution from the nonresonant graphs, and the 
solid lines are the total. 
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An interesting property of the quark loop amplitudes is the cancellation of the Higgs 
graphs with the box diagrams which can be related to the good high energy behavior of 
tt — > ZZ. The relative sign between the triangle and the box diagrams has been checked 
with the unitarity condition. We can cut the loop diagrams and obtain the tree processes 
tt — > ZZ and bb — > ZZ. Unitarity then requires a cancellation among tree diagrams at high 
energy. No such cancellation is required for the scalar loops between diagrams with and 
without Higgs bosons. 

Vs~ = 14 TeV, M A = 300 GeV 
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Fig. 6. Differential cross section of pp — > ZZ + X from gg — ► ZZ via both quark and 
squark loops, for rriq = 200 GeV, as a function of Mzz, for (a) tan/3 = 2 and (b) 
tan/3 = 10. Also presented is the destructive interference between diagrams with 
the H and the h for (c) tan (3 = 2 and (d) tan (3 = 10. 
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In Fig. 6 we have plotted the contributions to gg — > ZZ including both the resonant 
Higgs graphs and the nonresonant graphs for tan/3 = 2 and 10 including both the quark 
loops and the squark loops. The background from squark loops can be safely neglected 
around the Higgs pole region. 

The H resonance exhibits destructive interference just below the peak. This arises from 
interference with the contribution from the diagrams involving the light Higgs boson h. This 
feature is present for squark loops alone and also for quark loops alone, and should be present 
in nonsupersymmetric extensions of the Standard Model with extended Higgs sectors. In the 
MSSM it seems unlikely that this interference dip can be seen experimentally given the size 
of the underlying background. We note that the contributions to gg — > ZZ from squarks 
loops without Higgs bosons are also valid in nonminimal supersymmetric models since they 
do not involve contributions from the Higgs sector. 

V. HIGGS SEARCH AT THE LHC 

The basic features that arise from the numerical calculations are that the squark produc- 
tion cross section peak near the real production thresholds, \fs ~ 2m^. This enhancement 
should be smoothed somewhat after the inclusion of QCD corrections. The overall level of 
the squark loop contribution is much smaller than that of the quark loops, and even after 
including the interference between the two contributions we find that the squark loops can 
always be safely neglected in an overall estimate of the background level. 

In Fig. 7, we present the invariant mass distribution of ZZ including the Higgs signal 
from gluon fusion and the irreducible background from qq annihilation and gg fusion via 
quark and squark loops in pp collisions, with a rapidity cut \yz\ < 2.5, for nig = 200 GeV, 
tan/3 = 2 and tan/? = 10. For the smaller values of tan/? near 2, pronounced peaks appear 
at Mzz = Mff- For tan/? > 10, the Higgs peaks become almost invisible. 

The Higgs signal in the "gold-plated" mode H — > ZZ — > A£ is viable provided there are no 
supersymmetric decays of the Higgs boson JJ. The "silver-plated" mode H — > ZZ — > £ + £~z/z7 
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might also be used, and has a rate six times that of the gold-plated mode, but reducible Z-jet 
backgrounds might be a problem f2"3|,|3|j. If real squarks can be produced the Higgs width 
becomes larger, the Higgs peak in the invariant mass of ZZ is reduced and the event rate is 
too small to discover the Higgs. Even if we choose squarks with a mass such that the squark 
contribution peaks near the Higgs pole, the contribution to the continuum background from 
the squark loops is negligible. 

Vs~ = 14 TeV, |y(Z)| ^ 2.5 
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Fig. 7. Differential cross section of the Higgs signal in pp collisions with a rapidity 
cut \y z < 2.5|, from quark and squark loops, for rriq = 200 GeV, as a function of 
Mzz, for tan/? = 2 and 10; and, the irreducible background da/dMzz via qq and 
gg fusion. 

We expect that the squark loop background is also negligible for Higgs bosons with 
masses Mh,h < 2Mz, so that at least one of the Z bosons is virtual {Z*). In fact for the 
lighter supersymmetric Higgs, h, this is an attractive channel for discovery at the LHC [|J. 



The calculation gg — > ZZ* — > A£ has been performed for the quark loops in Ref. [24|, and 



that background falls rapidly for invariant masses decreasing below the ZZ threshold. We 
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expect the squark loops to be negligible in this case as well, but our calculation is applicable 
only for real Z bosons. 



VI. CONCLUSIONS 



We have calculated the contribution from squark loops to the process gg — > ZZ at the 
LHC. We find that they are substantially smaller than the contribution from quark loops, 
and can usually be safely neglected in computation of the continuum background in the 
Higgs search. In models with two CP-even Higgs bosons, we have also shown the amplitudes 
involving the heavy Higgs boson interferes destructively just below the resonance with the 
contributions from the light Higgs bosons. The expressions for the helicity amplitudes given 
in the appendix can be easily adapted to give the squark loop contributions to the processes 
gg — > 77, 77 — > 77 and 77 — > ZZ. The first process is relevant for Higgs searches at the 
LHC, while the last two are relevant for the 77 options of future linear e + e~ colliders. 
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APPENDIX 

In this appendix we give the formulas for the helicity amplitudes. We define helicity 
vectors for the process (12 — > 34) as 
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2s A Y 



m z (t — u) ~2 + t + m z — A Tn^-jj + m z \- 1 + m z 



-{t^u,(3 s ^ -A) 



m?F(t — u) 



-C (s) 



2s 4 

2 r (t- M )(t 2 + M 2 + 2F) 



-C(s) 



-B(s) 



4s 4 
s 2 s 2 



S(t 2 + M 2 )~ 



y 



(t - «) + A 



4y 

s(i — u + A s ) 



t 2 + u 2 - 2m% 
s 4 



2s 4 

sY(t — u — sp a ) 

2s 4 tiMi 



(26) 



A 



++++ 



771; 



D(s : t)2mi [mi + —*- 



s 4 



+c(0 



SS 4 ti 



s 2 (*i(*- «) +y) -2m 2 z ti 



+A -2m 2 



'(l + /3 s >-2m|' 
. 2s 4 " , 



gg^y 

S 2 S 4 



? ti(t-u) + y tiy 

S 4 ti ss 4 



,7771 



- 2- - 1 + (i <-> u) 



+£>(*, «)^— [(2S777? + y) (2SS 4 777 2 + S 2 y) + f3 s sY (y + 4S777 2 ) 



+C (s)4 



2s 2 s 4 

777|m| 

s 4 



18 



+ 



1 



m%{2Y - ss 4 ) + (3 s sY 



S4t\U\ 

where = s — 4m|, s 2 = s — 2m|, t\ = t — m|, U\ = u — m|, (3 S 



(27) 



s 4 /s, Y = tu — m|, 



and A = y—2m 2 z /sY. The squark mass-squared is denoted by Trig, where this is the 
eigenvalue (s) of the squark mass-squared matrix. Different masses for the left- and right- 
handed squarks can be accounted for by separating the overall factors in Eq. |l5|,and using 
two squark masses m~ and m~ in the expressions for the helicity amplitudes above. The 



factors B, C and D are the usual scalar integrals that result from the reduction [|Tl]?0 of 
the tensor integrals in the original Feynman diagrams. Expressions for these in terms of 
logarithms and dilogarithms can be found in Ref. ||. The other helicity amplitudes are 
related by the same substitutions required for the quark loops M 



A\ 1 A2A3A4 


— A* 


-A2— A3— A4 ' 




A++—(I3) 




-+(-/*) , 




A+++-(J3) 


= A++. 






A + — (P) 


= Af-H 


-+(/?) , 




A+-+({3) 


= Af-H 


-W) , 




{(5) = A ++0+ ({3) 


= A++. 


. (-P) = A++o-{ 


-p) 


= -A+- +(-P) 


= A+— 


-o(-P) = -A + -o- 


-((3) 



(28) 

The helicity amplitudes for top squark loop diagrams containing the Higgs bosons for 
the process gg — > ZZ are: 

2m\ sin a 



A Higgs 
• /i ++00 



a<,a 



sin 2 9 W cos 2 9 

S-2 



mz (a , > 
+ — cos(p + a) 



x^U + 2m 2 C (,)) 



raw sin (3 2 cos 9. 

cos(/3 — a) 



cos 6L 



s — m H + iTnTTiH 



+ 



2m? cos a 



m z 



mw sin (3 2 cos 9. 
x^(l + 2m 2 C ( S )) 



sin(/3 + a 



j cos9 w 
sin(P — a) 



2ml 



s - m 2 h + iT h m h 



(29) 
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A 



Higgs 



a, a 



sin 2 9,„ cos 2 9„ 



2m 2 sin a 



+ 



m z 



x(l + 2m 2 -C (s)) 



sin {3 2 cos ft,, 
cos(/3 — a) 



cos(/3 + a) 



s — mjj + iT H rriH 



+ 



2m 2 cos a; 



mw sin /3 2 cos ft, 
x(l + 2m 2 -C (s))- 



■ sin(/3 + a) 
sin(/3 — a) 



cos ft. 



s - m| + iT h m h 
and the contributions from the bottom squark are: 



A Higgs 
^++00 



a s a 



sin 2 6 W cos 2 ft. 



2m 2 cos a 



mw cos /3 2 cos ft. 



■ cos(/3 + a) 



+ 



s 2 „ o^y / \ \ cos(/5 - a) 
x^(l + 2m 2 C ( S ))- 

z 

2m 2 sin a 



s — + iT H mH 



+ 



m z 



mw cos (3 2 cos 0„ 
«2 „ , „ 2/ 



^Higgs 



2m 2 z 



a<,a 



;i + 2m-C ( S )) 



sin(/3 + a 
sin(/9 — a 



cos ft. 



s - m? + iT h m h 



l ++H 



sin 2 ft,, cos 2 9,, 



"h 
m z 



2m 2 cos a 



mw cos /3 2 cos ^ 
cos(/5 — a) 



cos(/3 + a) 



x(l + 2m 2 C ( S ))- 2 

s — m# + tL H^H 



+ 



2m 2 sin a 



+ 



m z 



mw cos /3 2 cos 0. 



sin(/3 + a) 



m z 
cos ft. 



x(l + 2m 2 C (s))- 



sin(/3 — a) 



cos 6 1 ,, 



cos ft. 



cos 9„ 



(30) 



(31) 



(32) 



s-m 2 h + iT h m h _ 

We remark here that the equations above involving transverse Z bosons can be reduced 
to the helicity amplitudes for gg — > 77 from squark loops not involving the Higgs graphs 
(the process gg — > h — > 77 has been investigated including squark mixing in Refs. f5|,|10|1). 
Also the contribution to photon-photon scattering 77 — > 77 can be obtained from that for 
gg — > 77 by changing overall factors. 
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